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Abstract —The problem of resolving the fine details of a signal 
from its coarse scale measurements or, as it is commonly referred 
to in the literature, the super-resolution problem arises naturally 
in engineering and physics in a variety of settings. We suggest a 
unified convex optimization approach for super-resolution. The 
key is the construction of an interpolating polynomial in the 
measurements space based on localized kernels. We also show 
that the localized kernels act as the connecting thread to another 
wide-spread problem of stream of pulses. 

I. Introduction 

Super-resolution is the task of estimating the fine details of a 
signal from its coarse scale measurements (see for instance Q], 
[ 0 , |3|). This problem arises in many physical and engineering 
problems as sensing systems have a physical constraint on the 
highest resolution the system can achieve. The aim of this 
paper is to suggest a unified convex optimization approach 
for super-resolution problems, based on the construction of 
interpolating polynomials and the existence of well-localized 
kernels. 

For simplicity, we consider signals of the form 

x{t) = y ^c m 5 trn , teM, (1.1) 

m 

where St is a Dirac measure, M C d > 1, is a compact 
manifold and T := {t m } is the signal’s support. This model 
can be extended to piece-wise polynomials given additional 
information (e.g. boundary condition, signal’s average). 

The information we have on the signal is its low spectral 
coefficients, namely its ’low-resolution’ measurements. In this 
paper, we consider three concrete examples. Suppose M is 
the ci-dimensional torus. In this case, we assume that the 
measurements are given as 

y = F n x , 

where Fjy projects x onto the space of trigonometric polyno¬ 
mials of degree N. That is, we have access solely to the low 
2 N + 1 low Fourier coefficients of x. This model reflects the 
fact that many measurement devices have limited resolution 
(e.g. diffraction limits in microscopes) . This problem can 
be solved using parametric methods such as MUSIC, matrix 
pencil and ESPRIT 0, 0, (6), Q. However, the robustness 
of these methods is not well understood. We mention that 
asymptotic results regarding the robustness of these methods 


were derived in (8), (9| and significant steps towards a non- 
asymptotic behavior have been taken recently ED, CD. 

Another closely-related example is the projection of x onto 
the space of algebraic polynomials of degree N. In this case 
we choose M as the manifold [—1, l] d . This problem arises, 
for instance, in spectral methods to numerically solve partial 
differential equations (see e.g. ca, on). 

The third example is low-resolution measurements of sig¬ 
nals lying on the bivariate sphere § 2 . In this case, we assume 
that a Dirac ensemble of the form (1.1 1 lies on § 2 , and we 
have access only to its low N spherical harmonics coefficients, 
which is the extension of Fourier analysis for signals on the 
sphere. This problem arises for instance in medical imaging 

oa, ca. 

The rest of the paper is organized as follows. The following 
section describes the convex optimization framework and its 
main ingredient - the construction of interpolating polynomi¬ 
als. Section |nl] reveals the intimate relations between localized 
kernels and super-resolution and demonstrates it on the three 
examples. Later on. Section IV is devoted to the problem 
of recovery from stream of pulses and its connection to the 
super-resolution problem. Section [V] shows several numerical 
experiments for super-resolution on the sphere. Ultimately, 
in Section [VI] we draw some conclusions and suggest future 
extensions. 


II. Convex Optimization Approach to 
Super-resolution 


In this paper we focus on a convex optimization approach 
for resolving signals from their low resolution measurements. 
We use the Total-Variation (TV) norm as a sparse-promoting 
regularization. In essence, the TV norm is the generalization 
of /' | norm to the real line (for rigorous definition, see for 
instance Section 1 in El and El). For signals of the form 
(U}, we have ||x|| T v = l c ™l- 

The main pillar of this framework is the following theorem 

El: 


Theorem II.1. Let x(t) = Y{ m c mSt m ,c m £ K, where T := 
{t m } C M, and M is a compact manifold in R d . Let 
be a linear space of continuous functions of dimension D in 
M. For any basis {Pk}k Li of Hd, let yk = ( x,Pk) for all 
1 < k < D. If for any signed sequence set {u m } £ {—1,1} 





there exists q £ YYjj such that 


q(tm) = U m , Mtm £ T , (H.l) 

\q(t)\<l,Vt£M\T, (11.2) 


then x is the unique real Borel measure satisfying 


min ||iH t’v subject to yk = (x,Pk), Y < k < D, 
xeM(M) 

(H.3) 

where A 4(M) is the space of signed Borel measures on M. 
By taking II£> as the space spanned by the eigenfunctions 


of the low-resolution measurement operator. Theorem II. 1 


implies that the super-resolution problem can be reduced to 
a polynomial construction problem. We mention that sim¬ 
ilar techniques, frequently referred to as dual certificate , 
are widely used in many sparse recovery problems, see for 
instance lfl9l . (20l . l2fl . 

From the algorithmic perspective, it is interesting to notice 
that in all three cases the infinite-dimensional TV minimization 
problem ( |II.3| > can be solved accurately (i.e. to any desired 
resolution) by algorithms with finite complexity which de¬ 
pends only on N. The recovery is performed by a three-stage 
algorithm involving the solution of a semi-definite program 
based on the dual problem of ( |II.3| > (can be solved using off- 
the-shelf software), followed by root-finding and least square 
fitting l22l . fl6ll . 11231 . Alternatively, one can discretize the 
problem and solve standard t\ minimization which converges 
to the solution of ( |lL3l ) as the discretization becomes finer 
[|24l . As we discuss later, these solutions are robust to noisy 
measurements under a separation condition. 


III. The Construction of Interpolating 
Polynomials 

As aforementioned, super-resolution problem can be re¬ 
duced to the construction of an interpolating polynomial which 
lies in the space spanned by the eigenfunctions of the low- 
resolution measurement operator. The existence of such a 
polynomial relies on two interrelated pillars. The first is the 
separation condition defined as follows. 

Definition III.l. A set of points T C M is said to satisfy the 
minimal separation condition with respect to the metric df, •) 

if 

A := min d(ti,tfi>p/N, 

tifyeTjijHj 


Consider a signal of the form of 0 defined on the 
circle [0,1] and its projection onto the space of trigonometric 
polynomials of degree N. In the time domain, the projection 
is given as 

y(t) = (x * D N )(t), 

where Djy(t) = e ^ kt * s the Dirichlet kernel of degree 

N. 

The key for super-resolving the signal is the existence 
the Fejer kernel Dn, a well-localized smooth super-position 
of Dirichlet kernels. Equipped with the wrap-around metric 
d(ti,tj ) = ||C — tj 11oo, the authors of lf22l showed that 
under the separation condition of Definition |III.1| there exist 
coefficients {a m } and {6 m } so that the polynomial 


q(t) = ^2 a mD N (t - t m ) + bmD^ (t - t m ), 

m 


satisfies ( |II. 1 ) > and ( |II.2| >. Hence, the TV minimization ( |II.3| ) re¬ 
covers x exactly in the univariate case. They also showed that 
similar polynomials can be constructed for higher dimensional 
signals. In this manner, the existence of the Fejer kernel is 
crucial for the ability to resolve a signal on the d-dimensional 
torus under the separation condition. In consecutive papers, 
the existence of the interpolating polynomial also serves as 
the basis for proving that the TV minimization results in a 
robust and localized solution in the presence of noise |25l , 

ESI. 


A similar result holds for the case in which the measure¬ 
ments are the projection of a signal of the form (1.1 1 on [—1,1] 
onto any basis spanning the space of algebraic polynomials of 
degree N. As was shown in lITSll . by considering the metric 
d{ti,tj ) = |arccos (fi) — arccos (fj)! one can construct the 
appropriate algebraic polynomial q{t) under the separation 
condition of Definition |III.1| and hence the recovery through 
TV minimization is guaranteed. Observe that in this case, the 
minimal required separation is space-dependent and reduced 
near the edges to the order of O (IV -3 / 2 ) . A consecutive paper 
showed that the recovery is also robust to noisy measurements 
E3- These results hold for the bivariate case as well. 

The same phenomenon occurs for signals on the bivariate 
sphere § 2 . Recall that spherical harmonics is the natural 
extension of Fourier analysis to signals on the sphere, and 
consider a Dirac ensemble on § 2 


®(0 = Y, Cm5 ^ {64 e a c s 2 . (III.l) 

m 


where v > 0 is a constant which does not depend on N. 

Along this work we argue that the separation condition 
is a sufficient condition. However, we emphasize that few 
works prove that the separation is also necessary, and without 
minimal separation the recovery can not be stable in the 
presence of superfluous noise E2, E3, ED- 

The second pillar is the existence of well-localized kernels 
in the space spanned by the eigenfunction of the low-resolution 
measurement operator. To explain this statement, we will 
demonstrate it through our three prime examples. 


Let Y n k{i), 0 < n < N, —n < k < n be an orthonormal 
basis of Vn (§ 2 ), the space of spherical harmonics of degree 
N. We assume that the information we have on the signal is 
its projection onto Vn (§ 2 ) 

2/n,fc = (x, Y Utk ), 0 < n < N, —n < k <n. (III.2) 

In the space domain, the projection onto Vn (§ 2 ) can be 
written as a spherical convolution, 

y(0=[ x(rj)K N {£-r))dS 2 (r)), 

J s 2 








IV. Super-resolution and Stream of Pulses 


where by the addition formula |[28ll 

N 

K n(^'V) = ^2 Y n,k(0 Y n,k(v) = ~^~ Pn > 3 ^'^), 

k=—N 


where Y n tk is the conjugate of Y n<k , and Pat ,3 (at) is the 
univariate ultraspherical Gegenbauer polynomial of order 3 
and degree TV. 

In l29l . it was shown that a smooth super-position of the 
kernel results in a well-localized kernel in the space 

of spherical harmonics of degree TV. We denote the localized 
kernel as _Fjv(£- 77 ), and by , £= 1,2 the partial rotational 
derivatives at £. Leveraging the localization of F N , it is known 
ED that there exist coefficients {a m }, {/3 m } and { 7 m} so that 
a polynomial of the form 


q(£) = V a m F N ((, • £ 


m 

2 £Tv(£, 


fulfils the requirements of Theorem II. 1 under the separation 
condition with the natural distance on the sphere d (£«,£/) = 
arccos (£j ■ £/■). Accordingly, we present the following theo¬ 
rem: 


Theorem III.2. m Let 5 = {£ m } be the support of a signed 
measure x of the form (III. 1). Let {Y n t k}n—o b e an y spherical 
harmonics basis for VnJW) and let y n ^ = (x , V n ,fe)> 0 < 
n < N, —n < k < n. If 2 satisfies the separation condition 
of Definition III.l with the metric d(£i,£j) = arccos (fii • £j), 
then x is the unique solution of 


min ||x|| T v subject to y n ,k = {x,Y nk ) 

xGA^(S 2 ) 

0 < n < TV, —n < k < n, 


(III.3) 


where ,M(S 2 ) is the space of signed Borel measures on S 2 . 

It has been shown in f23l (for a discrete version of ( |III.3| l) 
that the recovery error in the presence of noise is proportional 
to the noise level. 

We conclude this section with an interesting observation 
regrading non-negative signals (i.e. c m > 0). In this case, 
a sufficient condition for signal recovery is the existence of 
interpolating polynomial as in Theorem |II.1| but the constraint 
in (II.1 1 is replaced by a weaker constraint of q(t m ) = 1 for 
all t m £ T (see Theorem 5.1 in lfl 8 l ). Consider the case of the 
projection of signals on [ 0 , 1 ] onto the space of trigonometric 
polynomials of degree TV. In this case for any s < TV a 
polynomial of the form 


S 

q{t) = 1 - 2 ~ (s+1) (1 - cos (t - tm)) , 

m =1 


is a trigonometric polynomial of degree N, q(t m ) = 1 for all 
t m £ T and |g(f)| < 1 otherwise. Consequently, a clustered 
TV-sparse signal with non-negative coefficient can be recovered 
exactly. The same observation is noted for the other two cases 
with similar constructions (see fT 8 i l. 


In the previous section we stressed that the existence of 
well-localized kernels is crucial for resolving signals. The 
localized kernels bind the super-resolution problem to the 
problem of stream of pulses, namely recovery of the delays 
T := {t m } and weights {c m } from the measurements 

y(t) — ^ ( C-mdS-fj (t tm) 3 C m , t £ M, 

m 

where K(t) is a pulse shape (kernel) and I\ a {t) \= K{t/a) for 
some scaling parameter a > 0. An alternative representation 
of the problem is as y(t) = (x * K a ) (t) where 

ttt{f) — ^ ) C-mdtjn ; Y .— {fm} 3 
m 

and St is a Dirac measure. 

In a manner similar to Theorem |II. 1| the recovery problem 
can be reduced to construction of a special interpolating 
function, as follows: 

Theorem IV.l. I&0\l Let x{t) = c mSt m , c m G K where 
T := {t m } C R, and let y{t) = J R iV(f — s)dx ( s) for a L 
times differentiable kernel I\(t). If for any set {u m } £ {— 1,1} 
there exists a function of the form 

q(t)= [ ^ { %-t)dMs), 

Jr t _ 0 

for some measures {pi {t)}i =0 , satisfying 

q(t m ) = ttm 3 £ X) 

\q(t)\ < 1 , ft £ R\T, 


then x is the unique real Borel measure solving 
min I lx 11 tv subject to 

5GA4(R) 

r av.i) 

y(t) = / K(t — s)dx(s). 

Jr 

We note that here too, as in the super-resolution problems, 
the existence of q(t) relies on separation and localization. If 
the support T satisfies the separation condition with the metric 
d(ti, tj) = | ti — tj | and a = 1/TV, there exist coefficients {a m } 
and {b m } so that 


q{t) = ^2 a mKa ( t - tm) + b m K^ [t - t m ) ■ 


satisfies the requirements of Theorem IV. 1| and thus enabling 
perfect recovery. This holds if the kernel K{t) is localized. 
More precisely, the kernel should satisfy the definition of 
admissible kernel as follows: 


Definition IV.2. A kernel K is admissible if it has the 
following properties: 

1) K £ C 3 (R), is real and even. 

2) Global property: There exist constants Ce > 0,1 = 
0,1, 2,3 such that |( t) | < Ci/ (l + t 2 ). 

3) Local property: There exist constants e, (3 > 0 such that 
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Fig. 1: The mean recoery error (in log scale) as a function of v over 
20 simulations. To be clear, by error we merely mean the distance 
on the sphere between the true and the estimated support. 


a) K(t) > 0 for all |i| < e and K(t ) < K(e) for all 
1*1 > e > 

b) I\ ^ (t) < —fi for all |t| < e. 

Here too, the existence of the interpolating function guaran¬ 
tees a robust and localized recovery ESI, ED- Interestingly, 
in the non-negative case (i.e. c m > 0) the separation is 
unnecessary. In the presence of noise, the recovery error is 
proportional to the noise level, and depends on the number of 
spikes within any resolution cell of size vo (32l . 


V. Numerical Experiments on the Sphere 


To demonstrate our results we chose to use the less familiar 
example where we consider the recovery of a Dirac ensemble 
on the sphere (III. 1 1 from its low-resolution measurements 
( ]III.2| ). In this case, the TV minimization ( |I1I.3[ > can be solved 
to any desired accuracy by a three-stages algorithm consists 
of semi-definite optimization problem with O (A -4 ) variables, 
followed by root-finding and least square fitting. In case that 
the measurements are contaminated with bounded noise, it has 
been shown that the recovery error is proportional to the noise 
level (in the discrete setting) l23l . 

The following experiments were conducted in Matlab using 
CVX ES, which is the standard modelling system for convex 
optimization. The signals were generated in the following two 
stages: 

• Random locations on the sphere were drawn uniformly, 
sequentially added to the signal’s support, while main¬ 


taining the separation condition of Definition III. 1 


Once the support was determined, the amplitudes were 
drawn randomly from an iid normal distribution with 
standard deviation of SD = 10. 


Figure [T| presents a numerical estimation of the separation 
constant v (see Definition |HL1) . As can be seen, a separation 
constant of 27 t seems to be sufficient in the noise-free setting. 
This separation coincides with the spatial resolution of the 
projection of x onto Vjy El- 

Figure [2] presents an example for super-resolution on the 
sphere. The signal in Figure [2a] is the projection of the signal 
onto Do, whereas Figure[2b]presents the recovered signal. The 


N 




(a) The low resolution measurements. 



(b) The recovered signal /. 

Fig. 2: Super-resolution on the sphere for N = 10. The signal is 
presented on a grid for visualization only. 


N 

5 

8 

10 

Average error 

8.1267 x 10~ 5 

8.1826 x lCT 5 

9.0404 x 10 -5 

Max error 

2.163 x 10“ 4 

1.9 x 10” a 

3.3 x lO' 15 


TABLE I: The localization error for N = 5, 8,10. For each value of 
N, 10 experiments were conducted. To be clear, by error we merely 
mean the distance on the sphere between the true and the estimated 
support. 


maximal and average recovery errors for several values of N 
are presented in Table [I] 

In the noisy setting, we added an iid noise with normal dis¬ 
tribution and zero mean. Figure [3] presents the mean recovery 
error as a function of the noise standard deviation. We note 
that the error increases moderately with the standard deviation. 

VI. Discussion 

In this paper, we have presented a general framework for 
resolving robustly signals in various settings and geometries 
using convex optimization. Focalized kernels play a crucial 
role in the process. An important question is how general is 
this framework. That is to say, under which settings we can 
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Fig. 3: 10 experiments were conducted with N = 8 for each value 
of standard deviation. The figure presents the average recovery error 
as a function of the noise standard deviation. 


expect to resolve signal robustly using convex optimization. 
This topic is under ongoing research. 

We also showed that the localization principle relates super¬ 
resolution to other problems such as the recovery from stream 
of pulses. Recently, new results on the localization of the 
solution of ( |IV.1[ ) have been developed and tested on real 
ultrasound data Oil . It is of a great interest to find more 
applications where similar techniques could be applied (for 
recent related works, see for instance ED, ED). 

From the algorithmic perspective, the aforementioned super¬ 
resolution problems can be recast as finite dimensional prob¬ 
lems which involve the solution of a semi-definite program, 
root-finding and least square fitting. It will be interesting to 
look for the relations of these algorithms with ’traditional’ 
parametric methods such as MUSIC which also rely on root- 
finding. 

References 

[1] H. Greenspan, “Super-resolution in medical imaging,” The Computer 
Journal, vol. 52, no. 1, pp. 43-63, 2009. 

[2] S. C. Park, M. K. Park, and M. G. Kang, “Super-resolution image 
reconstruction: a technical overview,” Signal Processing Magazine, 
IEEE, vol. 20, no. 3, pp. 21-36, 2003. 

[3] M. Elad and A. Feuer, “Restoration of a single superresolution image 
from several blurred, noisy, and undersampled measured images,” Image 
Processing, IEEE Transactions on, vol. 6, no. 12, pp. 1646-1658, 1997. 

[4] P. Stoica and R. L. Moses, Spectral analysis of signals. Pearson/Prentice 
Hall Upper Saddle River, NJ, 2005. 

[5] Y. Hua and T. Sarkar, “Matrix pencil method for estimating parame¬ 
ters of exponentially damped/undamped sinusoids in noise,” Acoustics, 
Speech and Signal Processing, IEEE Transactions on, vol. 38, no. 5, pp. 
814-824, 1990. 

[6] R. Roy and T. Kailath, “Esprit-estimation of signal parameters via rota¬ 
tional invariance techniques,” Acoustics, Speech and Signal Processing, 
IEEE Transactions on, vol. 37, no. 7, pp. 984-995, 1989. 

[7] R. Schmidt, “Multiple emitter location and signal parameter estimation,” 
Antennas and Propagation, IEEE Transactions on, vol. 34, no. 3, pp. 
276-280, 1986. 

[8] H. Clergeot, S. Tressens, and A. Ouamri, “Performance of high res¬ 
olution frequencies estimation methods compared to the cramer-rao 
bounds,” Acoustics, Speech and Signal Processing, IEEE Transactions 
on, vol. 37, no. 11, pp. 1703-1720, 1989. 

[9] P. Stoica and T. Soderstrom, “Statistical analysis of music and subspace 
rotation estimates of sinusoidal frequencies,” Signal Processing, IEEE 
Transactions on, vol. 39, no. 8, pp. 1836-1847, 1991. 

[10] W. Liao and A. Fannjiang, “Music for single-snapshot spectral estima¬ 
tion: Stability and super-resolution,” arXiv preprint arXiv:1404.1484, 
2014. 


[11] A. Moitra, “The threshold for super-resolution via extremal functions,” 
arXiv preprint arXiv:1408.1681, 2014. 

[12] C. Canuto, M. Hussaini, A. Quarteroni, and T. Zang, Spectral Methods 
in Fluid Dynamics (Scientific Computation). Springer-Verlag, New 
York-Heidelberg-Berlin, 1987. 

[13] J. Shen, “Efficient spectral-galerkin method i. direct solvers of second- 
and fourth-order equations using legendre polynomials,” SIAM Journal 
on Scientific Computing, vol. 15, no. 6, pp. 1489-1505, 1994. 

[14] J. Tournier, F. Calamante, D. G. Gadian, A. Connelly et al., “Direct 
estimation of the fiber orientation density function from diffusion- 
weighted mri data using spherical deconvolution,” Neuroimage, vol. 23, 
no. 3, pp. 1176-1185, 2004. 

[15] S. Deslauriers-Gauthier and P. Marziliano, “Spherical finite rate of 
innovation theory for the recovery of fiber orientations,” in Engineering 
in Medicine and Biology Society (EMBC), 2012 Annual International 
Conference of the IEEE. IEEE, 2012, pp. 2294-2297. 

[16] T. Bendory, S. Dekel, and A. Feuer, “Exact recovery of non-uniform 
splines from the projection onto spaces of algebraic polynomials,” 
Journal of Approximation Theory, vol. 182, no. 0, pp. 7 - 17, 2014. 

[17] W. Rudin, Real and complex analysis (3rd). New York: McGraw-Hill 
Inc, 1986. 

[18] T. Bendory, S. Dekel, and A. Feuer, “Exact recovery of dirac ensembles 
from the projection onto spaces of spherical harmonics,” Constructive 
Approximation, pp. 1-25, 2013. 

[19] E. Candes, “Mathematics of sparsity (and a few other things),” 1CM 
2014 Proceedings, to appear, 2014. 

[20] G. Tang, B. N. Bhaskar, P. Shah, and B. Recht, “Compressive sensing 
off the grid,” in Communication, Control, and Computing (Allerton), 
2012 50th Annual Allerton Conference on. IEEE, 2012, pp. 778-785. 

[21] Y. De Castro and F. Gamboa, “Exact reconstruction using beurling min¬ 
imal extrapolation,” Journal of Mathematical Analysis and applications, 
vol. 395, no. 1, pp. 336-354, 2012. 

[22] E. J. Candes and C. Fernandez-Granda, “Towards a mathematical theory 
of super-resolution,” Communications on Pure and Applied Mathemat¬ 
ics, 2013. 

[23] T. Bendory, S. Dekel, and A. Feuer, “Super-resolution on the sphere 
using convex optimization,” Signal Processing, IEEE Transactions on, 
vol. 63, no. 9, pp. 2253-2262, May 2015. 

[24] G. Tang, B. N. Bhaskar, and B. Recht, “Sparse recovery over continuous 
dictionaries-just discretize,” in Signals, Systems and Computers, 2013 
Asilomar Conference on. IEEE, 2013, pp. 1043-1047. 

[25] E. J. Candes and C. Fernandez-Granda, “Super-resolution from noisy 
data,” Journal of Fourier Analysis and Applications, vol. 19, no. 6, pp. 
1229-1254, 2013. 

[26] C. Fernandez-Granda, “Support detection in super-resolution,” arXiv 
preprint arXiv:1302.3921, 2013. 

[27] Y. De Castro and G. Mijoule, “Non-uniform spline recovery from small 
degree polynomial approximation,” arXiv preprint arXiv:1402.5662, 
2014. 

[28] K. Atkinson and W. Han, Spherical harmonics and approximations on 
the unit sphere: an introduction. Springer, 2012, vol. 2044. 

[29] F. Narcowich, R Petrushev, and J. Ward, “Decomposition of besov and 
triebel-lizorkin spaces on the sphere,” Journal of Functional Analysis, 
vol. 238, no. 2, pp. 530-564, 2006. 

[30] T. Bendory, S. Dekel, and A. Feuer, “Robust recovery of stream of pulses 
using convex optimization,” submitted. 

[31] T. Bendory, A. Bar-Zion, S. Dekel, A. Feuer, and D. Adam, “Robust 
and localized recovery of stream of pulses with application to ultrasound 
imaging,” in preparation. 

[32] T. Bendory, “Robust recovery of positive stream of pulses,” submitted. 

[33] M. Grant and S. Boyd, “CVX: Matlab software for disciplined convex 
programming, version 2.1,” Mar. 2014. 

[34] B. Rafaely, “Plane-wave decomposition of the sound field on a sphere by 
spherical convolution,” The Journal of the Acoustical Society of America, 
vol. 116, no. 4, pp. 2149-2157, 2004. 

[35] R. Heckel, V. I. Morgenshtern, and M. Soltanolkotabi, “Super-resolution 
radar,” arXiv preprint arXiv:1411.6272, 2014. 

[36] C. Aubel, D. Stotz, and H. Bolcskei, “Super-resolution from short-time 
fourier transform measurements,” arXiv preprint arXiv:1403.2239, 2014. 





